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Fermilab fixed target experiment E791 obtained a sample of 1172 ± 61 events of D + — > iv~n + iv + and 848 ±44 
events of D+ -> n~TY + n + decays. We find respectively B(D + -> n~ tt+tt 4 ") / B(D + K~ty + ty + ) = 0.0311 ± 
0.0018+° °g26 and B ( D f -> vr + 7r-7r + )/B(L> + -» <^7r+) = 0.245 ± 0.028+° oi2 • Using a coherent amplitude analysis 
to fit the Dalitz plot of the D + — > 7r~7r + 7r + decay, we find strong evidence for a scalar resonance of mass 
478+23 ^ 17 MeV/c 2 and width 324+4Q ± 21 MeV/c 2 , compatible with what is expected for the isoscalar meson 
a . The D + — > <j(500)7t + fraction accounts for approximately half of all three-charged-pion decays of the D + 
. From the Dalitz plot analysis of the Df — > 7r~7r + 7r + decay events, we find significant contributions from the 
channels p°(770)7r+, p°(1450)tt + , /o(980)tt + , / 2 (1270)tt + , and / (1370)tt + . We also present new measurement of 
the masses and widths of the isoscalar resonances /o(980) and /o(1370). 



1. Introduction 

E791 data were produced by 500 GeV/c ir~ in- 
teractions in five thin foils (one platinum, four di- 
amond) separated by gaps of 1.34 to 1.39 cm. The 
detector, the data set, the reconstruction, and the 
resulting vertex resolutions have been described 
previously[ [l|. After reconstruction, events with 
evidence of well-separated production (primary) 
and decay (secondary) vertices were retained for 
further analysis. From the 3-prong secondary ver- 
tex candidates, we select a 7r~7r + 7r + sample with 
invariant mass ranging from 1.7 to 2.1 GeV/c 2 . 
For this analysis all charged particles are taken 
to be pions; i.e., no direct use is made of particle 
identification. 



2. D+ and D; 
ing Ratio 



relative Branch- 



We fit the spectrum of Figure [l] as the sum of 
D + and Df signals plus background. We model 
the background as the sum of four components: 
a general combinatorial background, the reflec- 
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tion of the D + — * K~ir + ir + decay, reflections of 
D° — > K~ir + plus one extra track (mostly from 
the primary vertex), and _D+ — > rj'ir + followed by 
i -> / o°(770)7, p°(770) -> 7T+7T-. We use Monte 
Carlo (MC) simulations to determine the shape of 
each identified charm background in the t7~tt + tt + 
spectrum [ || . We assume that the combinatorial 
background falls exponentially with mass. The 
fit finds 1172 ± 61 D+ events and 848 ± 44 D+ 
events. 

The branching fraction for Df — ► 7r~7r + 7r + rel- 
ative to that for Df — > 6tt + is measured to be: 



B(D+ -> TT--K+TT+) 



= 0.245 ± 0.028] 



(1) 



B(D+ -» 07T+ 

The branching fraction for D + — > 7r~7r + 7r + rel- 
ative to that of D + — ► K ~ir + ir + is measured to 
be: 

X + + Zk^C) = 0.0311±0.0018ta.(2) 

The first error is statistical and the second is sys- 
tematic. 

3. Dalitz plot formalism 

To study the resonant structure of these decays 
we consider the 1686 events with invariant mass 
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Figure 1. The tt~tt + tt + effective mass spectrum. 
The dotted line represents the D° — > K~ir + plus 
Df — > ?/7r + and the dashed line is the total back- 
ground. Events in the hatched areas at the D s 
and D + mass are used for the D s and D + Dalitz 
plot analyses. 

between 1.85 and 1.89 GeV, for the D + analysis 
and the 937 events with invariant mass between 
1.95 and 1.99 GeV/c 2 for the Df . Fig. 2a,b 
shows the Dalitz plot for these events. The hor- 
izontal and vertical axes are the squares of the 
7r + 7r~ invariant masses, and the plot has been 
symmetrized with respect to the two 7r + 's. 
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Figure 2. (a) The Df 7r"7r+7r+ Dalitz plot 
and (b) the Z? + — » 7r~7r + 7r + Dalitz plot. Since 
there are two identical particles, the plots are 
symmetrized. 

We fit each distribution shown in Figure |2| to 
a signal probability distribution function (PDF), 
which is a coherent sum of amplitudes corre- 
sponding to the non-resonant decay plus five 
different resonant channels, and a background 



PDF of known shape and magnitude. The 
resonant channels we include in the fit are 
p°(770)7r+, /o(980)tt+, / 2 (1270)tt+, / (1370)tt+, 
and p°(1450)7r+. 

We assume the non-resonant amplitude to be 
uniform across the Dalitz plot. Each resonant 
amplitude, except that for the /o(980), is param- 
eterized as a product of form factors, a relativistic 
Breit-Wigner function, and an angular momen- 
tum amplitude which depends on the spin of the 
resonance [ [| . 

For the /o(980)7r + we use a coupled-channel 
Breit-Wigner function, following the parameteri- 
zation of the WA76 Collaboration [ 01, 



BW 



1 



/o(980) 



ml v - ml + iTO (r 7r + T K ) 



with 

r w = 

and 



h^/mln/A-ml 



(3) 



(4) 



p 9K 
1 K = "7T 



y/ml„/A - m 2 K+ + - m 2 K ^j .(5) 



We multiply each amplitude by a complex co- 



efficient, Ci 



a s e s i. 



The fit parameters are the 



magnitudes, a,-, and the phases, <L, which accom- 
modate the final state interactions. The reported 
parameters values are obtained using the 
maximum-likelihood method. 

4. Df — > 7r~7r + 7r + Dalitz plot results 

The parameters of the /o(980) state, g n , gx, 
and mo, as well as the mass and width of the 
/o(1370), are determined directly from the Df — > 
7r _ 7r + 7r + decay events, floating them as free pa- 
rameters in the fit. The measured values Too = 
977±3±2MeV/c 2 , ffjr = 0.09 ± 0.01 ± 0.01 and 
g K = 0.02 ± 0.04 ± 0.03. For the / (1370) we find 
mo =1434 ± 18 ± 9 MeV/c 2 and r = 173 ± 32 ± 6 
MeV/c 2 . In both results, the first reported error 
is statistical and the second is systematic. The 
other resonance masses and widths are taken from 
the PDG[|. 

Table 1 shows the phases (6j) determined from 
the fit, and corresponding fraction for each decay 
mode. We calculate the decay fraction for each 
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mode 


relative phase 


fraction(%) 


/o(980)tt+ 


0° (fixed) 


57. ± 4. ± 5 


non-resonant 


(181. ± 94. ± 51.)° 


.5 ± 1. ± 2. 


p u (770)vr+ 


(109. ± 24 ± 5.)° 


6. ± 2. ± 4 


/ 2 (1270)tt + 


(133. ± 13. ± 28.)° 


20. ± 3. ± 1. 


/o(1370)7r+ 


(198. ± 19. ± 27.)° 


32. ± 8 ± 2. 


p u (1450)7r+ 


(162. ± 26. ± 17.)° 


4. ± 2. ± .2 



Table 1. Results of the Df Dalitz plot fits (sys- «> 
tematic error follows the statistical). 30 



amplitude as its intensity, integrated over the 
Dalitz plot, divided by the integrated intensity 
of the signal's coherently summed amplitudes. 
The first reported error is statistical and the sec- 
ond is systematic, the latter being dominated 
by the uncertainties in the resonance parame- 
ters, in the background parameterization, and 
in the acceptance correction. The /o(980)vr + is 
the dominant component, accounting for nearly 
half of the Df — > ir~ir + TT + decay width, followed 
by the /o(1370)vt + and /2(1270)vt + components. 
The contribution of / 9°(770)vt+ and ( o°(1450)7r+ 
components corresponds to about 10% of the 
yT~7r + 7r + width. We have not found a statisti- 
cally significant non-resonant component. The 
two m 2 + _ projections are nearly independent 
and the sum of them is shown in Fig. |^. 

To assess the quality of our fit absolutely, and 
to compare it with other possible fits, we de- 
veloped a fast-MC algorithm that simulates the 
£)+ — + 7T~vt + ?t + Dalitz plot from a given signal 
distribution, background, detector resolution and 
acceptance. For any given set of input parameters 
we calculated a \ 2 using the procedure presented 
in Ref . [ |j| . From x 2 and the number of degrees 
of freedom (y), we calculate a confidence level 
assuming a Gaussian distribution in \ 2 l v - The 
confidence level for the agreement of the projec- 
tion of our result onto the Dalitz plot with the 
data of table I is 35%. 

5. D + — > ?t 7t + ?t + Dalitz plot analysis re- 
sults 

In a first model for this channel, which we 
will refer to as Fit 1, the signal PDF includes 
the same resonances states used in the £>+ — > 
vt~vt + ?t + analysis. In the case of the /o(980) and 
/o(1370), we used the parameters obtained in the 
D + — » 7r~7r + 7r + analysis. In this model, the non- 
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Figure 3. projections for data (dots) and 

fast-MC (solid). The hashed area is the back- 
ground distribution. 

resonant, the p (1450)vT + and the p°(770)vT + am- 
plitudes dominate [ |^| . The qualitative features 
of this fit are similar to those reported by E691[ 
|andE687[|. 

We produce the x 2 distribution for the differ- 
ence in densities and observe a concentration of 
X 2 in the low 7r + ?r~ mass 2 (m 2 + _) region. The 
X 2 summed over all bins is 254 for 162 degrees of 
freedom (y), which corresponds to a confidence 
level less than 10 -5 , assuming Gaussian errors. 
We display the sum of the two m 2 „. in Fig. |4^i for 
the data and for the fast-MC. 

The low confidence level casts doubt on the va- 
lidity of the model used. While the projection 
of the MC onto the / k + tt~ mass 2 axis describes 
the data in the /9°(770) and /o(980) regions well, 
there is a discrepancy at lower mass, suggesting 
the possibility of another amplitude. To investi- 
gate the possibility that another -k + tt~ resonance 
contributes an amplitude to the D + — » 7r~vT + 7r + 
decay, we add a sixth resonant amplitude to the 
signal PDF. We allow its mass and width to float 
and assume a scalar angular distribution. This fit 
(Fit 2) converges and finds values of 478^ 2 g ± 17 
MeV/c 2 for the mass and 324^ 2 ± 21 MeV/c 2 
for the width, the first error is statistical and the 
second systematic. We will refer to this possible 
state as the <r(500). 

In Fit 2, the cr(500) amplitude produces the 
largest decay fraction, 46%, with a relatively 
small statistical error, 9%. The non-resonant 
fraction, which at (39 ± 10)% was the largest in 
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Figure 4. m„ w projection for data (error bars) 
and fast-MC (solid) . The shaded area is the back- 
ground distribution, (a) solution with the Fit 1, 
(b) solution with Fit 2. 

the original fit, is now only (8 ± 6)%. When 
we project this model onto the Dalitz plot, the 
X 2 /j/ becomes 138/162. The projection of this 
model onto the ir + ir~ invariant mass squared dis- 
tribution, shown in Fig. ^b, describes the data 
well, including the accumulation of events near 
0.2 GeV 2 /c 4 . 



mode 


relative phase 


fraction(%) 




(207. ± 8. ± 5.)° 


46. ±9. ±2. 


P u (770)tt+ 


0° (fixed) 


34. ± 3. ± 2. 


non-resonant 


(57. ± 20. ± 6.)° 


8. ± 6. ± 3. 


/o(980)tt+ 


(165. ± 11. ± 3.)° 


6. ± 1. ± .4 


/2(1270)tt+ 


(57. ± 8. ± 3.)° 


19. ± 3. ± 0.4 


/o(1370)7r+ 


(105. ± 18. ± 1.)° 


2.5 ± 1.5 ± 1. 


P u (1450)tt + 


(319. ± 29. ± 11.)° 


1. ± 1. ± .3 



Table 2. Results of the D + Dalitz plot fit, includ- 
ing the scalar er(500) resonance, (the systematic 
error follows the statistical). 

To better understand our data, we also fit it 
with vector, tensor, and toy models for the sixth 



(sigma) amplitude, allowing the masses, widths, 
and relative amplitudes to float freely. The vector 
and tensor models test the angular distribution of 
the signal. The toy model tests the phase vari- 
ation expected of a Breit-Wigner amplitude by 
substituting a constant relative phase. These al- 
ternative explanations of the data fail to describe 
it as well the scalar with a regular Breit Wigner 
option 

6. Conclusions 

Using a Dalitz plot analysis of the three pion 
decay of the Df, we find a little but signifi- 
cant contributions from the channels p°(770)ir + , 
p°(1450)7r + . We measure the mass and width of 
the /o(980) and we have not found a statistically 
significant gx component in the width. Using 
the same analysis for the D + mass region, we 
find strong evidence that a scalar resonance with 
mass 478 1 24 ± 17 MeV/c 2 and width 324±^ ± 21 
McV/c 2 produces a decay fraction « 50%. These 
results are discussed in references [ ||, ||, || . 
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